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ELECTRON  DRIFT  IN  A  LINEAR  MAGNETIC  WIGGLER  WITH  AN 
AXIAL  GUIDE  FIELD 


Recently,  there  has  been  much  interest  in  free  electron  lasers  (FEL's) 
which  use  spatially  varying  magnetic  fields  to  modulate  a  relativistic 
electron  beam  (REB).*-^  Three  types  of  magnetic  wigglers  have  been  widely 
used:  helical  (usually  produced  by  a  bifilar  helical  current  winding**^), 

Q 

linear  (produced  by  an  array  of  permanent  magnets  or  by  linearly  alternating 
windings^),  and  radial*®  (produced  by  a  series  of  spaced  conducting  or 
ferromagnetic  washers**  immersed  in  an  axial  field  or  by  a  series  of 
alternating  coils  ) .  Of  these,  only  the  first  two  produce  a  perpendicular 
field  on  axis  and  are  therefore  suitable  for  use  with  a  small  diameter  beam  of 
solid  cross-section.  In  this  paper,  we  will  analyze  the  propagation  of  a 
solid,  high  current  REB  through  linear  and  helical  wigglers  with  a 
superimposed  axial  guide  field. 

FEL  experiments  fall  into  two  categories  depending  on  the  beam  current. 
Compton  regime  FEL's*’^1^  have  used  high  energy  (10's  of  MeV) ,  low  current 
(<  1  A)  beams  while  Raman  FEL's  have  used  lcwer  energy  (~  1  MeV),  high 
current  (>  1  kA)  beams.  Helical  wigglers  have  been  used  in  both  current 
regimes,  but  until  recently^  linear  wigglers  were  used  only  in  the  Compton 
regime.  In  this  case,  the  linear  wigglers  consist  of  permanent  magnets  and 
there  is  no  guide  field.  However,  for  high  current  beams,  a  guide  field  is 


required  to  contain  the  beam.  Although  there  are  advantages  to  a  linear 


wiggler  from  the  standpoint  of  ease  of  assembly  and  versatility  (e.g., 
changing  the  periodicity  or  tapering  the  period  and/or  field  amplitude  for 

i  o 

efficiency  enhancement  )  we  will  show  that  there  is  no  equilibrium  for  off- 
axis  particles  when  a  beam  is  propagated  through  a  linear  wiggler  in  an  axial 
field.  The  particles  drift  out  of  the  wiggler  unless  additional  focusing 


forces  are  provided. 
Maniucript  submitted  April  2, 1982. 
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Various  authors  have  considered  particular  cases  of  electron  motion 


through  magnetic  wigglers.  These  have  typically  been  single  particle  (low 

current)  calculations.  Blewett  and  Chasman^  considered  motion  of  high  energy 

electrons  (~  24  MeV)  through  a  helical  wiggler  and  found  stable  helical  orbits 

with  superimposed  betatron  oscillations.  Friedland  has  treated  the  case  of 

lower  energy  electrons  (  ~  300  keV)  in  an  idealized  radially  uniform  helical 

wiggler  with  a  superimposed  axial  guide  field.  He  showed  that  various  stable 

trajectories  were  possible  and  derived  stability  criteria  relating  the  allowed 

wiggler  and  guide  field  strengths  to  beam  energy  and  wiggler  period  X  . 

w 

These  "stable"  regions  are  given  by 
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where  BQ  is  the  axial  guide  field,  »  eB^fm  is  the  corresponding  cyclotron 
frequency,  v  is  the  electron  velocity,  and  k  ■  2tt/Xw.  Physically,  these 
conditions  stem  from  the  resonance  in  the  perpendicular  velocity  of  an 
electron  in  the  wiggler  and  guide  fields:^ 
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Freund  and  Drobot^  have  considered  this  case  further  and  also  conclude  that 
stable  trajectories  with  nearly  constant  axial  velocities  and  relatively  large 


wiggler  amplitudes  are  possible  when  <<  ck.  This  is  consistent  with 
condition  (2a). 


The  present  analysis  employs  a  relatively  simple  computer  code  which 
solves  the  equations  of  motion  of  an  electron  in  any  electric  and  magnetic 
field  configuration  using  a  fourth-order  Runge-Kutta  method.  Self  electric 
and  magnetic  fields  are  calculated  by  assuming  that  the  electron  is  at  the 
edge  of  an  azimuthally  symmetric  beam  of  current  I,  so  that  the  self  fields 
can  be  written  as  (mks) 


"r  2ire  rv 
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where  r(t)  is  the  electron  radius.  The  axial  self  fields  are  neglected. 


For  the  particular  cases  considered  here,  the  external  magnetic  field 
consists  of  a  solenoidal  field  produced  by  a  15.3-cm-I.D.  x  2-m-long  solenoid 
together  with  a  wiggler  that  begins  in  the  uniform  portion  of  the  solenoidal 
field.  The  wiggler  field,  which  may  be  either  helical  or  linear,  rises 
adiabatically  over  ten  periods  and  then  oscillates  with  constant  amplitude. 
The  envelope  enclosing  the  wiggler  amplitude  is  given  by 
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Here,  z  is  Che  axial  distance  from  the  beginning  of  the  wiggler.  This 
variation  closely  fits  the  linear  wiggler  used  in  the  NRL  induction  linac 
FEL.  The  field  on  axis  from  this  wiggler  is  plotted  in  Fig.  1  along  with  the 
envelope  equation  (5). 

The  linear  wiggler  field  components  are 

B  -  b(z)B  cosh  kx  cos  kz 
X  w 

By  -  0  (6) 

J l 

B  *  -b(z)B  sinh  kx  sin  kz* 
z  w 

where  Bw  is  the  peak  wiggler  field  on  axis.  For  the  NRL  linear  wiggler,^ 

which  has  -  3  cm  and  winding  layers  spaced  by  3.2  cm,  these  expressions  are 

valid  for  r  <  1  cm. 

««** 

We  will  be  primarily  concerned  with  relatively  low  energy  (~1  MeV) ,  high 
current  (~kA)  beams.  Of  particular  interest  are  the  effects  of  self-fields, 
guide  field  and  wiggler  amplitude,  and  initial  beam  conditions  corresponding 
to  a  field- immersed  (v0o  *  0)  or  a  shielded  source  (PQq  *  0),  where  v0q  is 
the  lntitlal  azimuthal  electron  velocity  and  P0q  the  initial  canonical 
angular  momentum.  If  the  axial  self  magnetic  field  is  neglected,  P0q  »  0 
implies  vg  ■  erBz/(2ym)  in  a  uniform  field  Bz. 

First,  we  consider  a  beam  with  y  m  2.2,  I  -  250  A,  BQ  -  2  kG,  B,,  *  1  kG, 
and  Aw  -  3  cm.  Note  that  Eq.  (2a)  would  require  Bw  <  3  kG  for  stability  in 
this  case  with  a  helical  wiggler.  The  x-y  trajectories  of  electrons  injected 
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with  v^q  •  0  and  rQ  ■  0.4  cm  into  a  linear  wiggler  are  shown  in  Fig.  2.  Each 
particle  initally  begins  to  E  x  I  drift  in  the  self  fields,  but  when  it 
reaches  a  region  of  large  wiggler  field  it  begins  drifting  in  the  y- 
direction.  In  each  case  the  electron  reaches  an  imaginary  wall,  located  at  r 
*  1.2  cm,  after  only  partially  traversing  the  wiggler.  It  should  be  noted 
that  the  assumptions  used  to  calculate  the  self  fields  become  invalid  as  the 
beam  distorts.  Consequently,  to  determine  if  these  self  fields  are 
responsible  for  this  drift,  the  calculations  are  repeated  with  I  »  0  (Fig. 

3) .  A  third  trajectory  is  also  plotted  for  Bw  *  0.5  kG.  The  electron  that 
originates  on  the  y-axis  is  well  confined,  but  electrons  off  the  y-axis  again 
drift  to  the  wall,  with  a  velocity  much  higher  for  Bw  **  1  kG  (<v^>  «  .047  c) 
than  for  ^  3  .5  kG  ^vd  ^  c)  • 

This  behavior  can  be  explained  by  the  increasing  gradient  in  the  wiggler 

field,  and  consequently  the  emergence  of  a  significant  axial  field  component, 

as  the  distance  from  the  y  axis  increases.  Although  the  present  drift  arises 

from  the  gradient  in  B„  and  is  in  the  direction  of  VB  x  B  ,  it  is 

w  wo 

quantitatively  different  from  the  usual  guiding  center  approximation  because 

the  field  variation  over  one  gyroperiod  is  so  extreme.  For  example,  when  kx  ■ 

% 

.8  (x  »  .4  cm  in  the  present  case) ,  Bz  varies  from  +  Bw  to  -  Bw  over  one 
period.  However,  the  physical  mechanism  is  the  same  as  for  the  usual  gradient 
drift;  i.e.,  the  gyroradius  in  the  part  of  the  orbit  where  Bz  is  a  minimum,  or 
where  1x1  is  a  minimum  in  Fig.  (3),  is  larger  than  where  Bz  is  a  maximum.  The 
addition  of  self  fields  merely  imposes  an  additional  Ex?  rotation  on  the 
obit,  so  that  an  electron  that  originates  on  the  y-axis  (and  is  therefore 

l 

confined  when  1*0)  begins  to  drift  into  a  region  of  increasing  B 

z 

Therefore,  it  actually  is  "lost"  sooner  than  an  electron  on  the  x-axis  which 
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>  >  J) 

initially  E  x  B  drifts  into  a  region  of  smaller  Bz 


An  approximate,  empirical  expression  for  the  drift  has  been  found  which  is 

in  quite  good  agreement  with  the  code  results  when  £2  <  Q  <  kv  (below  the 

wz  o  z 

cyclotron  resonance).  In  terms  of  the  wiggler  gradient,  the  expresson  is 


»  I-/  V(BS2  ^o 

Vd“2  “o  2(bVXbo 


(7) 


which  has  the  same  form  as  the  usual  grad  B  drift  but  is  quantitatively 

different.  The  single  particle  equations  of  motion  are  ■  Vy^z * 

v  «  -  v  £2  +  v  £2  ,  and  v  ■  v  SI  .  Assuming  v  *>  constant  and  that  the 
y  zxxz  z  yx  °z 

gyration  velocities  are  v  »  v.  sin  kz  and  v  -  v.  cos  kz ,  it  follows  that 
oj  y  ly  x  lx 
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i  i 

for  Bx  »  B^  .  These  assumptions  are  reasonably  valid  for  kx  <  .8  and  for 

BQ  sufficiently  far  from  resonance  and  Bw  small  enough  that  v^  «  vz  . 

Furthermore,  to  insure  that  cosh  kx  *  constant  over  the  orbit,  we  restrict 

v ,  «  v,  *  v .  (i.e .  ,  £2.  <  kv_)  .  Then  using  Eqs.  (8)  and  (6)  in  Eq.  (7)  , 

lx  ly  l  o  z 

we  obtain 
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Table  1  compares  drift  velocities  for  various  cases  from  the  code  with  1*0 
to  those  from  Eq.  (9).  In  general,  the  agreement  is  very  good. 


Clearly,  the  gradient  drift  is  more  severe  when  Bw  is  a  considerable 

fraction  of  BQ.  However,  if  Bw  is  held  constant  at  1  kG  and  Bq  increased  to  4 

kG,  v^(y  *  2.2)  is  not  significantly  reduced.  This  relative  insensitivity 
to  Bq  is  due  to  the  increased  v^  ,  and  hence  a  larger  v^,  as  BQ  approaches  the 

cyclotron  resonance  field  Br.  One  consequence  of  this  large  drift  near 

resonance  is  that  it  limits  the  degree  of  gain  enhancement  achievable  through 
the  magneto-resonance  effect.^’^ 

If  Bq  sufficiently  exceeds  Br,  the  drift  can  be  quite  small.  For 
y  *  2.2  and  *  3  cm,  the  resonant  magnetic  field  B^  *  7  kG.  With  BQ  *  8 
kG  and  Bw  *  1  kG,  the  particle  still  drifts  to  the  wall  as  shown  in  Fig.  4a. 
The  drift  is  now  in  the  opposite  direction  to  that  below  resonance 
since  v^  changes  sign  when  Bz  >  Br<  This  has  the  effect  of  changing  the  phase 
of  Vy  oscillations  with  respect  to  those  of  Bz;  i.e.  ,  v^  is  positive  (for  x  > 
0)  when  Bz  is  a  minimum,  so  that  the  particle  drifts  in  the  +y  direction. 

When  Bq  is  increased  to  10  kG  for  this  case,  the  drift  is  small  enough  that 
the  electron  remains  confined  for  >  30  periods  as  shown  in  Fig.  (4b).  The 
confinement  remains  very  good  when  the  self  fields  of  a  500  A  beam  are  added. 

In  principle,  operation  of  FEL  experiments  with  BQ  >  Bf  is  possible  and 
has  been  demonstrated  with  a  high  current,  y  *  3.5  beam  in  a  helical 
wiggler.^  However,  competing  processes  such  as  the  cyclotron  oaser^ 
interaction  can  produce  large  radiated  powers  at  frequencies  close  to  those  of 
the  FEL  interaction  in  this  beam  energy  and  magnetic  field  regime. 


10 


Fig.  4  -  Electron  trajectories  in  linear  wiggler  with  I  ■  0,  Y  “  2.2, 

Bw  *  1  kG,  and  rQ  *  .4  cm.  a)  BQ  -  8  kG,  b)  BQ  »■  10  kG. 
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Consequently,  analysis  of  experimental  results  Is  more  difficult.  Also, 

arbitrarily  large  guide  fields  are  not  possible  with  a  magnetically  shielded 

diode  simply  because  the  electrons  will  be  mirrored  by  the  field.  The 

equilibrium  beam  radius^  in  such  a  case  depends  only  on  y>  I,  BQ,  and  the 

beam  emittance  e,  so  for  particular  beam  parameters  suitable  values  of  BQ  are 

1  ft 

limited.  The  required  field  in  kG  for  a  matched  beam  radius  R  in  cm  can  be 
written 


B  - 


1.36  I 
R23Y 


11.56  e 


n 


(10) 


where  I  is  in  kA  and  is  the  normalized  emittance  in  rad  cm.  For  example, 
if  en  -  ,14  tt  rad-co  (about  the  lowest  value  expected  for  a  thermionic  cathode 
beam  with  I  =  750-1000  A)1®,  a  Y  *  2.2  beam  with  R  *  .3  (.5)  cm  requires 
BQ  **  5.5  (2)  kG  for  I  <  800  A.  In  this  regime,  the  beam  is  emittance 
dominated  so  that  the  required  field  is  relatively  insensitive  to  I.  In 
principle,  smaller  radius  beams  could  be  used  with  larger  BQ,  thereby  doubly 
reducing  v^.  However,  experimentally  this  is  very  difficult  at  high  current 

levels. 

To  analyze  the  effect  of  a  shielded  diode  on  propagation  through  the 

wiggler,  we  repeat  the  above  calculations  with  an  initial  Vg  corresponding 

to  P„  ■  0.  Note  that  I  ■  1.75  kA  in  this  case,  which  is  the  current  required 
wo 

for  constant  radius  propagation  in  only  the  solenoidal  field  with  these 
initial  conditions.  As  shown  in  Fig.  5,  the  electron  propagates  at  nearly 
constant  radius  until  the  wiggler  amplitude  becomes  large  enough  that  the 


gradient  drift  begins  to  dominate.  Then  the  particle  is  lost  just  as  in 

the  v.  ■  0  case. 

0o 

The  behavior  of  an  off-axis  electron  in  a  linear  wiggler  and  axial  guide 

1  3 

field  should  be  compared  to  that  when  a  helical  wiggler  is  used.  We 
approximate  the  helical  wiggler  field  by^ 

■  b(z)Bw(-sin  kz  [l  +  i  k2(3x2  +  y2)]  +  i  k2  x  y  cos  kz} 

(H) 

B^  ■  b(z)Bwfcos  kz  [l  +  ^  k2(x2  +  3y2)]  -  ^  x  y  sin  kz} 

B^  *  -kb(z)Bw[l  +  -g  k2(x2  +  y2)  ](x  cos  kz  +  y  sin  kz). 

This  expansion  of  the  true  Bessel  function  expression  for  Bh  is  valid 

for  kr  <  1.  Friedland's  treatment*'*  of  electron  propagation  in  this  case 

assumes  a  radially  uniform  wiggler,  and  his  stability  condition  (Eq.  2a)  is 

not  stringent  enough  when  the  radial  variation  is  included.  For  example, 

Friedland  finds  stable  orbits  when  y  ■  1.587,  X  *  4  cm,  v.  *  0,  r  -  0, 

w  io  o 

BQ  *  1.26  kG,  and  Bw  *  1.04  kG,  so  that  Eq.  2a  is  barely  satisfied,  and  we  can 
duplicate  his  results  if  we  remove  the  radial  variation  from  B*1.  However, 
with  the  B*1  given  in  Eq.  11,  we  find  that  the  wiggler  field  must  be  reduced  to 
~625  G  to  obtain  stable  orbits. 

Although  the  radial  dependence  in  B*1  does  narrow  the  allowable  range  of 
operating  parameters,  stable  orbits  with  <r>  *  cons't  are  achievable  with  a 

helical  wiggler  in  a  guide  field.  Electron  trajectories  in  a  helical  wiggler 
for  the  same  conditions  as  used  previously  with  a  linear  wiggler  are  shown  in 

15 


Figs.  6  and  7.  Figure  6  is  to  be  compared  to  Fig.  2  and  Fig.  7  to  Fig.  5.  In 
both  cases,  the  electron  is  well-confined.  It  is  interesting  to  note  that  in 
Fig.  6,  the  electron  born  at  (x,y)  ■  (.4,0)  initially  E  x  B  drifts  in  the  + 

9  direction,  but  then  reverses  direction  as  the  wiggler  amplitude 
increases.  Since  grad  B*1  is  radial,  the  grad  B*1  drift  is  in  the  -  9 
direction  and  does  not  lead  to  the  beam  expansion  observed  with  the  linear 
wiggler. 

The  linear  wiggler  drift  described  here  imposes  additional  constraints  on 
the  parameters  of  an  FEL  experiment.  Obviously,  the  beam  radius  nust  be  kept 
as  small  as  possible  to  minimize  particle  loss  from  the  edge  of  the  beam. 
Preliminary  experimental  results  by  our  group  with  a  field  Immersed,  apertured 
source  indicate  that  particle  losses  can  be  kept  acceptably  small  in  this 
way.  Also,  if  y  is  large  enough  that  Br  »  BQ  »  Bw>can  be  satisfied  for 
relatively  large  Bw,  then  the  drift  can  be  kept  small  while  achieving 
acceptably  large  vw. 

Finally,  it  should  be  noted  that  the  drift  arises  from  the  asymmetry  of 

the  linear  wiggler  and  the  corresponding  absence  of  focusing  forces  in  the 

direction  perpendicular  to  the  wiggler  field.  Therefore,  it  should  be 

possible  to  stabilize  the  drift  by  imposing  an  additional  focusing  force  in 

that  direction.  For  example,  preliminary  results  indicate  that  electron 

propagation  through  a  "square”  or  symmetrized  linear  wiggler  is  very  stable. 

£ 

Such  a  wiggler  has  an  additional  component  B^  -  cosh  ky  cos  kz  and  a 

# 

corresponding  addition  to  B  of  -Bu  sinh  ky  sin  kz.  An  electron  trajectory 
through  such  a  wiggler  is  shown  in  Fig.  8. 
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Electron  trajectories  in  symmetrized  linear  wiggler  with  I 

■  .4  cm,  and  v 


2-2,  B0 


2  kG,  B 


1  kG,  r 


0 


In  conclusion,  a  gradient  drift  has  been  shown  to  exist  for  a  linear 
wiggler  in  an  axial  guide  field.  The  drift  can  be  substantial  with  small  or 


large  beam  current  in  some  parameter  ranges,  for  a  wide  range  of  initial 
conditions.  However,  the  advantages  of  a  linear  wiggler  are  sufficient  in 
many  cases  to  either  limit  operation  to  a  "stable"  parameter  regime  or  to 
impose  additional  focusing  forces  to  stabilize  the  drift. 

The  authors  gratefully  acknowledge  useful  discussions  with  C.A. 
Kapetanakos,  H.  Freund,  and  C.M.  Tang. 
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